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ABSTRACT 

With the installation of the Cosmic Origins Spectrograph on the Hubble Space Tele¬ 
scope, measurements of the metal content of the low redshift intergalactic medium 
(IGM) are now available. Using a new grid-based model for diffuse gas coupled to the 
SAGE semi-analytic model of galaxy formation, we examine the impact of supernova 
feedback on the pollution of the IGM. We consider different assumptions for the reheat¬ 
ing and ejection of gas by supernovae and their dependence on galaxy circular velocity 
and gas surface density. Where metals are present, we find the most likely metallicity 
to be —1.5 <logiQ(Z/Z0)< —1.0 a.t z = 0, consistent with both observations and more 
sophisticated hydrodynamic simulations. Our model predicts that the regions of the 
IGM with the highest metallicities will be near galaxies with M* ~ Mq and 

in environments of densities ^ lOx the mean. We also find that 90% of IGM metals 
at z = 0 are ejected by galaxies with stellar masses less than 10 ^°'^^M©. 

Key words: 


1 INTRODUCTION 

Diffuse gas fills the Universe, surrounding and connecting 
galaxies, the densest peaks of the cosmic web. This filamen¬ 
tary gas accounts for over half of the baryons at z = 0 
| |Dave et al.| 1999)[Shull et al.||2012[ ), but is much less dense 
than the gas occupying galaxies and therefore more difficult 
to observe. Nevertheless, the structure of the intergalactic 
medium (IGM) has been observed in Ly-a emission ( |Can-| 
[talupo et al.| 2014 Martin et al. 20141, Ly-a absorption 
I Lynds|1971| etc.), and X-rays (Werner et al.|2008 Ma et al. 


galaxies has been performed with hydrodynamic simulations 


|2009[ [Planck Collaboration et al.|2013[ ). Future instruments, 
such as the Square Kilometer Array, will also be able to see 
the IGM in faint 21-cm emission ( Takeuchi et al.|2014 |. 

Cosmic gas is not necessarily pristine hydrogen and he¬ 
lium left over from the Big Bang. As galaxies evolve, they 
eject gas from their halos back into the diffuse regions of the 
cosmic web through various feedback processes, including 
supernovae and active galactic nuclei (AGN). The ejected 
gas includes metals created during the life cycle of stars, 
pushed into the IGM by outflows from galaxies. 

For example, using the Cosmic Origins Spectrograph 


on the Hubble Space Telescope, Shull et al. (20141 have mea¬ 
sured the metallicity of the IGM to be ~ O.IZq, a number 


in agreement with most theoretical estimates (see e.g. Cen 


et al.||200T| [Fang fc Bryan||2001| [Chen et al.||2003l |Oppen 


heimer et al.|2012[), usin^arge number of sight lines. 

Historically, modelling the cosmic gas in and in-between 


(e.g. Popping et al. 2009 Dave et al. 20101. Hydrosimula¬ 


tions are numerically detailed and accurate, yet computa¬ 
tionally costly. Their detail makes them ideal for studying 
individual systems or small volumes where one may wish 
to explore the complexities of galaxy formation. Baryons 
can also be traced outside of halos, allowing modellers to 
study the IGM as well. Their cost, however, makes it diffi¬ 
cult to run statistically significant numbers of simulations. 
This means they are prohibitive for exploring rare events, to 
sample unbiased distributions, or to demonstrate the statis¬ 
tics of a particular phenomenon. 

A complementary method is to run an N-body dark 
matter-only cosmological simulation first, then in post¬ 
processing couple an analytic model approximating the 


baryonic physics of galaxy formation (e.g. Croton et al. 


20061. This “semi-analytic” approach averages out many of 
the small-scale details, but produces integrated properties 
for galaxies and halos that can be directly compared against 
observed survey data. Semi-analytic models are computa¬ 
tionally inexpensive to run and, for an equivalent run-time, 
can produce orders-of-magnitude more systems to study. 
Furthermore, the increase in speed allows large areas of pa¬ 
rameter space to be explored. 

Traditional semi-analytics follow the bound portions of 
the N-body simulation, treating the cosmic web only as it 
falls onto the galaxy. To look at the IGM, several exten- 
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Figure 1. Diffuse dark matter particles in a 1.25 h~^ Mpc slice of 
milli-Millennium at 2 : = 0. The white circles are the halos, scaled 
to their respective virial radii. 


sions have been implemented using N-body simulations as 
a base. Most of these lack galaxy formation and are devel¬ 
oped to address a single problem, such as creat ing mock HI 
and ^Hell observations ([Takeuchi et al. 20141, finding the 


ionisation structure resulting from galaxies (e.g Kim et al. 
201^ [Park eTaTIpOldl, or measuring the light cone effects 


on the 21-cm signal ( Datta et al.|2014 ). 

We propose a generalized method for modelling the 
IGM that works in a similar fashion to traditional semi- 
analytic models. Rather than producing data focused only 
on a particular science question (e.g. ionisation, hydrogen 
density, etc.), our new model creates a full data set of gas 
in the intergalactic medium, which is coupled in a generic 
way with a semi-analytic galaxy formation model. We use 
the semi-analytic model SAGE (Groton et al. 2006 Cro¬ 
ton et al., in prep) to inject gas and metals into the IGM 
from galaxies through feedback, in addition to the primor¬ 
dial “unprocessed” gas already in the IGM that we infer 
from the unbound dark matter in our N-body simulation. 
This combined IGM+galaxy model, called SAGE With All 
Matter Included (SWAMI), bridges the gap between semi- 
analytic models of galaxy formation and their more expen¬ 
sive cousins, hydrosimulations. 

In this paper we introduce this new model and use it to 
consider the importance of supernova feedback on the metal 
content of the intergalactic medium. Using SWAMI, we test 
various scenarios of reheating and ejecting gas and metals 
from halos. 

This paper is organised as follows: in Section]^ we dis¬ 
cuss the different components of SWAMI and the results of 
our fiducial model. In Section|^ we extend the fiducial model 
to test different reheating and ejection scenarios. Einally in 
Section]^ we discuss our results, how they compare to hy¬ 
drosimulations and observations, and their implications for 
future observations. This is then summarised in Section]^ 


2 MODELLING GALAXIES AND THE 
DIFFUSE COSMIC GAS 

Modelling the diffuse cosmic gas requires three separate 
stages. The first part is the dark matter only N-Body 
simulation, for which we use the milli-Millennium Run, a 


small er version of the Millennium Simulation (Springel et al. 
2005[) developed for testing purpose^ and described in Sec¬ 


tion [2T] 

The second part is the semi-analytic model of galaxy 
evolution we use to create galaxies within the dark mat¬ 
ter halos of the N-Body simulation. For this, we use the 


Semi-Analytic Galaxy Evolution (SAGE) code (Croton et al. 
|2006[ Cr oton et al., in prep), discussed more fully in Sec- 
tion 12.21 

The third part is a new semi-analytic model of the dif¬ 
fuse gas outside of halos, called SAGE With All Matter 
Included (SWAMI), using the unbound particles from the 
N-Body simulation. We discuss this part in Section [2.3| 


2.1 milli-Milennium: The Dark Matter 


We use the smaller milli-Milennium rather than the full Mil¬ 
lennium simulation because the particle data is available for 
the former, in addition to the halo data needed for the galaxy 
evolution model. The milli-Millennium simulation was run 


using the WMAPl cosmological parameters (Spergel et al. 
20031 in a periodic box of side-length 62.5h~^ Mpc, with a 


mass-per-particle of 8.6 x 10°h~ Mq. 

We superimpose a (100)^ cell grid onto the box, giving 
each cell a side length of 625h~^ kpc and an average of 19.68 
particles (representing a total average mass per cell of 1.7 x 
Mq). Local density contrast (p/p) is calculated for 
each cell using the number of dark matter particles located 
within, and this is used as our definition of environment. 
For reference. Milky Way-type systems at « = 0 contained 
entirely within a cell have a density contrast of about 6, not 
including any unbound dark matter also in the cell. 

The particles are then divided into bound and un¬ 
bound with a friends-of-friends plus gravitational unbinding 
procedure using the halo finder SUBFIND ( [Springel et al.| 
20011. The bound particles are assigned to halos, which are 


then used in SAGE to make galaxies. The unbound par¬ 
ticles remain binned in the cells, which SWAMI uses to 
model the diffuse gas. An example of this is seen in Fig¬ 
ure]^ which shows a 1.25h~^ Mpc (2 cell) deep slice of the 
milli-Millennium simulation. The white circles indicate halos 
(bound particles), scaled to their respective virial radii. The 
background structure shows the unbound particles binned 
into 625h~^ kpc cells, which we assume traces the diffuse 
cosmic gas. 


2.2 SAGE: The Galaxy Model 


Semi-analytic models have evolved considerably since they 
were first proposed by White & Frenk {19911, although they 
still follow many of the same assumptions. Every halo hosts a 


^ Full specs of the Millennium simulation can be found 
at http: //gavo .mpa-garching .mpg. de/Millenniuiii/Help?pag6= 
simulation 
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Figure 2. (left) Number of galaxies as a function of stellar mass and environment. The solid (dashed) line is a histogram of the 
environment (stellar mass) distribution, (right) Average mass of gas in the IGM in a cell as a function of stellar mass and environment. 
Both figures assume the fiducial SAGE feedback model. 


galaxy and each galaxy is made up of several different bary- 
onic reservoirs, including stars, cold gas (in the disk), hot gas 
(in the halo), and even ejected gas that has been pushed out 
of the halo. At each timestep, material cycles through the 
reservoirs according to analytical prescriptions of infall, gas 
cooling, star formation, galaxy mergers, and feedback pro¬ 
cesses. For our galaxy formation model we use SAGE, which 
is an update of the popular Croton et al. (2006 1 model, ex¬ 
cept for the modifications described below. 

For the creation and movement of metals, SAGE follows 
the prescriptions outlined in De Lucia et al. (20041. Metals 
are created by supernovae and pushed into the disk at a rate 
proportional to the star formation rate, Am*: 


Quasars affect the high mass end of the stellar mass function, 
whereas supernovae have more of an impact on the low mass 
end. 

One of the main sources of uncertainty in semi-analytics 
(and understanding galaxy evolution in general) is the mech¬ 
anism of feedback and its dependency on the evolving 
galaxy. To this end, we consider several different models for 
reheating and ejecting gas from the disk and halo with su¬ 
pernova winds. Below we describe our fiducial model. Later 
in Section 13 we consider several different models for reheat¬ 
ing and ejecting gas from the galaxy and halo. These form 
a key part of our results with SWAMI. 


Amr 


= Y^-Am*. 


( 1 ) 


We assume a yield, Y^, of 2.5%. As in many semi-analytic 
models, metal yield is a parameter of the model rather 
than an IMF-based calculation. Our value corresponds to 
a Salpeter initial mass function with an upper mass limit 
of 6 OM 0 . In contrast, De Lucia et al. (2004) use a yield of 
4.5%, which is similar to a Ghabrier IMF over the same mass 
range ( Madau fc Dickinson|2014 1. 

The energy output and timescales of the supernova pa¬ 
rameters we use are based mostly on core collapse super¬ 
novae, which happen soon after star formation and are re¬ 
sponsible for most of the alpha elements that are produced. 
Type la supernovae, in contrast, contribute to the heavy 
element population and occur on a delayed timescale from 
the star formation. Since elliptical galaxies tend to be older 
and therefore have had more of a build up of metals from 
Type la supernovae, the assumption of a constant yield di¬ 
rectly proportional to the star formation rate does not hold 
as strongly. Our model predictions therefore are most robust 
with alpha elements produced in late type galaxies. 

Once created, the metals then cycle through the various 
reservoirs - stars, cold gas, hot halo gas, and ejected gas 
- in a similar manner to the gas. Gas and metals can be 
ejected from a galaxy’s halo by both supernova and quasars. 


2.2.1 The Fiducial Galaxy Feedback Model 

To set a baseline for our results we use the default SAGE 
model described in Groton et al. (in prep.), including all 
parameter choices. We refer the interested reader there for 
the more general details. Below we focus on the parts of the 
model most relevant to the topic of the current work, namely 
the modelling of supernova feedback and the movement of 
metals beyond the galactic disk into the halo and IGM. 

Feedback processes, including supernovae and AGN, are 
crucial to the galaxy formation process. As stars form, a 
certain percent of the mass is automatically recycled back 
into the disk. The birth of very short-lived stars quickly 
result in supernovae that reheat the surrounding cold gas 
and push it out of the disk into the hot halo reservoir. We call 
this “reheating”, and assume that the amount of gas that 
is reheated (Amreheat) is proportional to the star formation 
rate and the disk mass loading factor, Cdiak: 


Amreheat — Cdisk ' Am*. 


( 2 ) 


In our fiducial model, Cdisk is a constant value with no de¬ 
pendence on the galaxy which hosts the supernova. This 
simple characterisation comes from early observations that 


measure ediak~ 


1 - 


5 dMartinll 1999 1. 


The energy released during a supernova explosion is not 
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completely consumed by the reheating of gas from the disk. 
It also drives winds that can eject hot gas from the halo. We 
assume conservation of energy, where the supernova energy 
left over from reheating the gas from the disk is used to ki- 
netically drive gas from the halo into the IGM. Conservation 
of momentum has been ruled out as a likely mechanism for 
supernova ejection ( Murray et al.|[^005 K 

From conservation of energy, the amount of energy 
available to eject gas, Ameject, is 


where 


and 


AEeject — AEsn AEreheat, 


AEsn — -ehaioAm*VsN 


1 2 

AEreheat ~ ^^^teheatM vir ■ 


(3) 


(4) 


(5) 


Ehaio is the efficiency with which the supernova energy trans¬ 
forms into kinetically driven winds and Vsn is the charac¬ 
teristic velocity of the wind, taken to be 680 km s“^. We 
assume the reheated gas reaches virial equilibrium in the 
halo on a time-scale shorter than the time resolution of the 
simulation (<250-350 Myr), where the halo virial velocity is 
Vvir. Assuming the escape velocity is ~ Vvir, Equations j^- 
[ 5 ] combine to: 


* _ ( 0 /^ VsN V \ 

ZAnieject — I ^halo I "y . I ^disk I 


Am* 


( 6 ) 


We define the halo mass loading factor to be ehaio= 
£haio(^V^) ’ ''^hich simplifies our expression for Ameject to: 


Ameject — (chalo €disk) ' Am*. 


(7) 


As with the heating of the gas to the virial temperature, we 
assume ejection due to supernovae (and subsequent partial 
reincorporation) occurs within a timestep of the model. 

After running the fiducial SAGE model on the halo 
merger trees of the milli-Millennium simulation we produce 
a catalogue containing 26,292 galaxies more massive than 
10^h~^ Mq at z = 0. In the left panel of Figurewe show 
the distribution of these galaxies as a function of stellar mass 
and environment (local cell density contrast, as described in 


Section 2.1 1 Many galaxies fall below the stellar mass reso¬ 
lution limit of our model, which is taken to be M* ~ 

M 0 . Galaxies less massive than this usually have dark mat¬ 
ter halos only just above the halo mass resolution limit of 
the simulation and are unlikely to have fully developed. Fig¬ 
ure shows the expected trend of more massive galaxies 
living in denser environments, with an additional grouping 
of low mass galaxies, just above the resolution, in environ¬ 
ments about 100 X the mean density. The solid (dashed) line 
along the x-axis (y-axis) highlights the distribution of galax¬ 
ies as a function of environment (stellar mass). 


particles are traced with pristine gas and that all of the dif¬ 
fuse metals are in the immediate region of the galaxies that 
have ejected them. While this is in contrast to hydrosimula¬ 
tions which have revealed metals not only outside the halos 
of galaxies (e.g. Brook et al.|2012 Stinson et al.|20T2l, but 


also in relatively empty regions of space (e.g. Oppenheimer 
et al. 20121, the assumption that the gas that has not yet 


been processed through the halo is pristine is a reasonable 
approximation at low redshift. To address high redshift mix¬ 
ing, including a lower limit of the metallicity of infalling gas, 
a more complex model is needed, and will be the subject of 
future work. 

We start by overlaying an analytic mapping between 
cosmic baryons and the unbound particle grid of the sim¬ 
ulation. We assign hydrogen and helium to the unbound 
matter according to the prescription in [Madau (2002 I: 

( 8 ) 


nn = {pcrit/mii){l - y)flf,(l -|- zf 


where pcrit, ran, Y, and Qb are the critical density, mass of 
hydrogen, the primordial abundance of helium (0.24), and 
the baryonic density (0.04), respectively. This gives us a local 
hydrogen content of 


Mh = nH(a:) x v^eii = nn)! -k 5) x Vceii 
= 1.12 X 10“®cm"^nt,h^(l + (5) X Vceii 


(9) 


where 5 is calculated using the total density of the cell, and 
Vceii is the volume of the cell. The amount of diffuse gas is 
adjusted by subtracting off the mass of the galaxy. 

To partially couple the galaxy model to the diffuse 
model, we add the ejected reservoirs of gas for each fiducial 
SAGE galaxy to the IGM in the cell where each galaxy is lo¬ 
cated. The median amount of diffuse gas surrounding galax¬ 
ies as a function of stellar mass and environment is shown 
in the right panel of Figure As expected, the amount of 
IGM surrounding a galaxy is much more dependent on en¬ 
vironment than on stellar mass, with IGM increasing with 
increasing density. The amount of IGM in a cell reaches its 
maximum at densities between 30 and 300 times the mean. 
Cells with higher densities than this contain mostly bound 
particles and thus do not have much IGM. 

The ejected metals predicted by each SAGE galaxy and 
ejected due to supernova and AGN feedback are added to 
the cells in the same manner as the ejected gas. We assume 
perfect mixing within the cell but none outside. Previous 
studies have found that ejected matter is recycled through 
the halo and galaxy an average of three times over the life of 


the galaxy (Oppenheimer & Dave 20081, and gas typically 
takes ~ 1 Gyr to cycle through the galactic fountain (|Brook| 


et al. 20141, so our mixing model, so our mixing assumptions 


are a reasonable first-order approximation. Future work with 
a fully coupled model will be able to address the spread of 
metals out of the immediate region of each galaxy/halo sys¬ 
tem and follow the history of metals in the IGM as galaxies 
move between cells. 


2.3 SWAMI: The Diffuse Cosmic Gas Model 

The IGM in our diffuse gas model is considered in two parts 
- the baryons that trace the unbound dark matter particles 
from the simulation and the ejected gas reservoir of each 
galaxy modelled by SAGE. We assume that the unbound 


3 EXPLORING SUPERNOVA WINDS AND 
METAL EJECTION 

The presence of metals in the hot galactic halo modulates 
the rate at which gas cools onto the galaxy, which in turn 
affects the star formation rate. By changing the mechanism 
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Figure 3. The three panels show the nine combinations of heat¬ 
ing and ejection models, used to calculate the ratio between 
ejected matter (Ameject) and amount of new stars formed (Am*), 
as shown in Eqn.[^ For all three panels, Rl, R2, and R3 are the 
dark, medium, and light lines, respectively. El, E2, and E3 are 
the solid, dashed, and dotted lines in the top, middle, and bottom 
panels, respectively. 


by which metals and gas are injected into and ejected out 
of the halo, we alter not only the metal content of the IGM 
but also the galaxy life cycle itself. Thus, the model assumed 
for gas reheating and ejection will have consequences for the 
wider galaxy population, although the details remain ill un¬ 
derstood. It is therefore worth exploring a range of possibil¬ 
ities and using these to compare with existing hydrosimula¬ 
tions and observations. In this section we focus on reheating 
(from the disk) and ejection (from the halo) separately, us¬ 
ing modihcations to the fiducial model disk and halo mass 
loading factors. We then examine the consequences of these 
modifications for selected galaxy and IGM properties. 


3.1 Gas Reheating And Ejection Models 


To test the effect of different reheating assumptions on the 
metallicity of the IGM, we vary the disk mass loading factor, 
edisk, according to various models suggested in the literature. 
We consider three in this paper: (Rl) our fiducial model of a 
constant euisk from [Croton et ah] | |2006p, (R2) a d isk circular 
velocity dependent euisk as used by Dave et al. (2010|), and 


(R3) a cold gas surface density dependent Cdisk from [Lagos 

eTahlpoT^: 


Cdisk — 3.0 

_ 150 km s~ 

^disk 


^disk — 


(r 


6xlOlS>M0Mpc-2 




Rl 

R2 

R3 


( 10 ) 

Rl: The simplest reheating model, Rl (our fiducial 


model), is commonly used by in semi-analytics and is moti¬ 
vated by the early work of Martin (19991, who showed that 
edisk~ 1 — 5 for galaxies over a wide range of circular veloc¬ 
ities, albeit with a large scatter. 

R2: Reheating model R2 is borrowed from a common 
implementation of galactic winds in hydrosimulations, such 


as by Oppenheimer & Dave (20061 and Hirschmann et al. 


(20131. A more physically motivated scaling than Rl, re¬ 
cent observations have indicated the disk mass loading fac¬ 
tor scales with circular velocity. For example, [Bordoloi et al.[ 
(20141 find that the equivalent width of galactic outflows 
increases with stellar mass for galaxies of M*< 

Mq at z ~ 1. We use halo Vmax as a proxy for disk circular 
velocity in our modelling. 

R3: Reheating model R3 is derived from dynamical 
modelling of superbubble expansion driven by supernovae 
( [Lagos et al.[[20l^ . Unlike Rl and R2, R3 takes its input 
from the SAGE galaxy properties directly (namely cold gas 
mass and stellar mass) and is therefore the most sensitive 
to the model assumptions and tunings (although note that 
we keep the same fiducial tuning for all results presented in 
this paper). 


As seen in Equation the pollution of the IGM from 
outflows is dependent on both euisk and thaio. To that end, 
we also explore the halo mass loading factor, thaio- We test 
models that vary both the ejection efficiency thaioj ^nd the 
supernova velocity scaling. The three versions of thaio we use 
are: 


^ _ ^0 f 680 km s ^ ^ 0 _ r, o tpi 

€halo ^halo I J ’ ^halo 

Chaic = eLio ^ 0.3 E2 

680 km ^ ^0 _ n i 170 

) 5 ^halo — U.l edisk 

( 11 ) 

El: Ejection model El assumes a fixed value of Chaio) 
set at 0.3 (our fiducial model), which represents a constant 
efficiency of supernova energy being transformed into kinetic 
winds. It also assumes a constant supernova wind velocity. 

E2: Ejection model E2 enhances the El by adding a de¬ 
pendence on the circular velocity of the galaxy. Here, we as¬ 
sume the supernova winds are not constant from one galaxy 
to another, but instead affected by the host galaxy. This is a 
similar assumption to our 2nd reheating model and a varia¬ 
tion of the models used in e.g. [Oppenheimer fc Davi| ( [2006[ ), 
etc. 

E3: Ejection model E3 adds a dependence on the disk 
mass loading factor through 0.1 Cdisk- This model 

comes from the assumption that the mechanism by which 
the supernovae reheat the gas from the disk is the same as 
the mechanism by which the gas is ejected from the halo. 
Physically, this model implies a strong interaction between 
the galaxy and its halo and that the amount of gas ejected 
depends only on the amount of gas that is reheated. Since 
hydrosimulations do not specifically eject gas from halos, 
this is a semi-analytic approximation for giving a particle a 
kick from a supernova and letting it travel. 

In Figure we plot the behaviour of the nine reheat¬ 
ing/ejection configurations. The three panels show the ratio 
of the mass of ejected gas to the mass of new stars formed 


_ 0 
Chalo — 
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during a timestep (equivalent to the difference between Chaio 
and Cdisk, Equation]^, with top, middle, and bottom panels 
showing El, E2, and E3, respectively. For all figures in this 
paper, dark, medium, and light coloured lines correspond to 
Rl, R2, and R3, respectively. Likewise, solid, dashed, and 
dotted lines refer to El, E2, and E3, respectively. RlEl and 
R1E3 are identical, as the fiducial value of Chaio E O.lx the 
fiducial value of tdisk. Since R3 is dependent on the cold gas 
and not directly on Vmax or Vvir, the R3 lines are from an 
empirical fit to the fiducial model and are meant for com¬ 
parison only. 

W e use the same convention as [Oppenheimer fc Dave| 
120081 and restrict E2 to a maximum Vmax of 226 km s“^. 
This corresponds to Mvir= 3 x Mq and limits the 

amount of matter that can be ejected from large halos. We 
include the factor of 3 so that galaxies at this maximum have 
the same supernova velocity as the other ejection models. 
There is a sharp cutoff for Rl with all ejection models at 
circular velocities corresponding to Mvir— lO^^/i"^ M0. 


3.2 Stellar Mass Function 

Before examining the IGM, we start by quantifying the ef¬ 
fects of our different reheating and ejection models on the 
galaxy population using the galaxy stellar mass function 
(SMF) at a = 0. Note that we do not retune SAGE between 
runs, only switch around the feedback model combinations. 
This provides a clean comparison of the impact each has on 
the model galaxy properties. The results are shown in Fig¬ 
ure]^ which, like Figure]^ has the three ejection models in 
three different panels. The green lines are the different mod¬ 
els with dark, medium, and light green showing reheating 
models Rl, R2, and R3, respectively. The solid, dashed, and 
dotted lines show ejection models El, E2, and E3, respec- 
region is the observed z = 0 SMF from 


Differences between the model combinations can mostly 
be seen at the low mass end of the SMF. For a constant ejec¬ 
tion efficiency and supernova velocity, (El, left panel), there 
is very little difference between the reheating model combi¬ 
nations, with all three matching the observations. On the 
other hand, galaxies that assume a circular velocity scaled 
ejection model (E2, middle panel) significantly overproduce 
the low mass population for all reheating models, with R3E2 
sitting slightly closer to the observed SMF compared to Rl 
and R2. Ejection model E3 (far right panel) is the most de¬ 
pendent on the reheating model. Here, R1E3 is the same as 
RlEl and is hence a good fit to the data, while R2E3 slightly 
over-predicts the SMF at 10 ^°Mq and under-predicts it be¬ 
low. R3E3 sits below the observed SMF at all masses. 

None of the feedback model combinations have a signifi¬ 
cant effect on the high mass end of the stellar mass function. 
This is to be expected, as star formation at the high mass 
end is typically suppressed by AGN feedback, which we are 
not varying. The high mass end is slightly under-predicted 
for all models because we are using the milli-Millennium 
simulation, which does not have a significant number of high 
mass galaxies. 

The nine combinations that are shown in Figure can 
be directly related to the stellar mass functions in Figure 
Gas that has been ejected must first fall back onto the 
halo before it can re-cool and settle onto the disk and form 


tively. The shaded 


Baldry et al. (2008 


stars. Therefore, models that eject less gas at lower masses 
(e.g. E2, middle panel) will produce more galaxies in that 
same mass range. Likewise, R3E3 ejects the most gas at low 
masses, so it produces the fewest low mass galaxies. An em¬ 
pirical fit to the fiducial model has edisk(R3) oc Vmax”^'^"^, 
making it more than an order of magnitude larger than 
edisk(R2) (oc Vmax“^) for sub-L* galaxies. 

Interestingly, R3 produces measurably fewer low mass 
galaxies than the other models do when the supernova wind 
speed is dependent on the galaxy’s circular velocity (E2), 
but still more than observations show, despite ejecting the 
least amount of gas in the mass range of interest (light blue 
line in the middle panel of Figure]^. This is due to Cdisk 
being very large at low masses and the supernova velocity 
being dependent on the circular velocity of the galaxy (and 
therefore the total available energy being much lower than in 
El). Most of the energy is therefore expended on reheating 
the gas from the disk, leaving both less energy to eject the 
gas and less cold gas in the disk to form stars right away. 
The reheated gas does not make it out of the halo, but it still 
has to cool back onto the disk to make stars, suppressing the 
SMF relative to the other reheating models. 


3.3 Total Metals in the IGM 

3.3.1 ... vs. Stellar Mass 

Figurej^shows the total mass of the metals ejected into the 
IGM by all galaxies of a certain mass. As in Figure the 
dark, medium, and light coloured lines are reheating models 
Rl, R2, and R3, respectively. The solid, dashed, and dotted 
lines are ejection models El, E2, and E3, respectively. 

For all models, the net ejected metals peaks at roughly 
M* « Mq. For the fiducial ejection model, El, our 

fiducial and velocity-scaled reheating models, Rl and R2, 
peak at slightly above 10 ^°Mq, but the cold gas scaled 
reheating model, R3, is slightly below. Reheating model R3 
produces fewer metals in the IGM than the other reheating 
models with the same ejection, while Rl and R2 produce 
similar amounts of metals. R1E2 and R2E2 differ by a factor 
of 1.5, but are still considerably closer to each other than 
the are to R3E2. R1E3 is not shown because it is the same 
as RlEl. Despite these differences, for all models at z = 0, 
90% of the metals in the IGM have been ejected by galaxies 
of M*< 10 ^°'®®Mq. 

Sub-L* galaxies are responsible for most of the met¬ 
als ejected into the IGM, and the amount of metals pro¬ 
duced and ejected are proportional to the star formation 
rate (Equations and [^. It would, therefore, stand to rea¬ 
son that the models that produce the most sub-L* galaxies 
would eject the most metals into the IGM. For example, 
for our reheating dependent ejection model (E3), the super¬ 
bubble derived reheating method (R3) produces the least 
amount of stars (right panel. Figure]^ and ejects a consid¬ 
erably lower amount of metals than our fiducial model. 

We find a different situation with the circular velocity 
scaled supernova wind ejection model (E2). All three reheat¬ 
ing models produce an excess of sub-L* galaxies (middle 
panel, Figure]^, indicating a high Am*, but have both the 
lowest and highest amounts of ejected metals (dashed lines. 
Figure]^. The variation in metal ejection for E2 models is 
due to the other part of the ejection equation (Equation]^, 
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logioM* (Mg) logioM* (Mg) logioM* (Mg) 


Figure 4. The stellar mass functions of the nine models. Left to right, the panels are constant El (solid lines), E2 (dashed lines), and E3 
(dotted lines). R1 (dark green lines), R2 (medium green lines), and R3 (light green lines) are shown in each panel. The shaded regions 
are the observed z = 0 SMF from [Baldry et al.| l |20d8l l. This figure uses h = 0.72 to accurately compare observations and theoretical 
models which have different dependencies on h. 


Ehaio - Ediak- The middle panel of Figure shows that, for 
the mass range in question, Ehaio-Edisk for E2 is highest when 
coupled with R2 and lowest when coupled with R3. This is 
reflected in the scalings in Figure 

Similarly, our fixed halo ejection efficiency and fixed 
supernova velocity ejection model (El) produces the same 
amount of stars regardless of the reheating model (left panel, 
Figure]^, indicating that Am* is roughly constant, so any 
variation in the amount of ejected metals relies entirely on 
the Ehaio- Edisk- The top panel of Figure shows the super¬ 
bubble driven reheating model (R3) is slightly suppressed 
at sub-L* masses compared to the other reheating models, 
so the R3E1 model (light solid line, Figure]^ puts the least 
amount of metals into the IGM for that ejection model. 


3.3.2 ... vs. Environment 

Figure]^ shows the total metals in the IGM as a function of 
the environment. As described in Section [2Tj environment is 
the density contrast of the dark matter particles in the cell. 
Like the above figures, the dark, medium, and light coloured 
lines are reheating models Rl, R2, and R3, respectively. The 
solid, dashed, and dotted lines are ejection models El, E2, 
and E3, respectively. 

Our results show that the mass of metals ejected into 
the IGM peaks in environments between 10 and 30 x the 
mean density, which are the typical environments of small 
groups. The ejection model E2 tends to peak at the lower end 
of this environment range (~ 10 x p) for all three reheating 
models. El and E3 peak in the middle (~ 18 x p) for Rl and 
R2, and at a higher density for R3E3. 90% of ejected metals 
are found in environments less dense than 40 x p. 

Galaxy mass scales roughly with density (Figure [^, 
with larger galaxies evolving in denser environments. Nev¬ 
ertheless, environment is much more closely related to halo 
mass than stellar mass. Comparing models as a function of 
environment is akin to comparing them as a function of halo 
mass. The models show very little variation in their slopes 
at lower densities, especially compared to the scatter be¬ 
tween the models at low stellar masses in Figure There 



Figure 5. Total ejected metals as a function of stellar mass. The 
dark/medium/light pink lines are reheating models 1, 2, and 3, 
respectively. The solid/dashed/dotted lines are ejection models 1, 
2, and 3, respectively. 

is more variation between the models at high densities, but 
their shapes are remarkably similar, differing only in their 
scaling. 

Environment, and by extension halo mass, is more im¬ 
portant for finding where the ejected metals are (in regions 
between 10 and 40 x p). The distribution of densities where 
metals are found in the IGM, i.e. the shape of the distribu¬ 
tion, is not model dependent, and the bulk of IGM metals 
will be found near small groups. The biggest difference be¬ 
tween the models is in the net amount of ejected metals, 
and because of this scaling, a given environment can have a 
vastly different metal content, depending on the model. 

Total metal mass in the IGM is more dependent on 
ejection models than on reheating models. For all ejection 
models, Rl and R2 have similar total mass and the same 
peak density. R3 has a much lower total mass than Rl and 
R2 for a given ejection model. For a given reheating model. 
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Figure 6. Total ejected metals as a function of environment. The 
dark/medium/light pink lines are reheating models 1, 2, and 3, 
respectively. The solid/dashed/dotted lines are ejection models 1, 
2, and 3, respectively. 

the various ejection models are much less tightly grouped, 
with no strong pattern emerging. 

3.4 Ejected Gas Metallicity 

Figurej^shows the average metallicity of the IGM as a func¬ 
tion of stellar mass and environment for three of the models 
(RlEl, R2E2, and R3E3), showing bins with at least 5 galax¬ 
ies. This figure is analogous to Figure but is now coloured 
by average IGM metallicity in the bin rather than number 
of galaxies or average surrounding IGM gas mass. The av¬ 
erage metallicity peaks at slightly sub-L* galaxies (10.0 < 
logj^Q(M*/h“^ Mq) < 10.5) in environments ~ lOx the av¬ 
erage density, with R2E2 peaking at slightly lower densities 
and stellar masses. The metallicity distributions of the mod¬ 
els that are not shown closely resemble the distribution of 
the same ejection model. 

Overall, the average metallicity in a stellar 
mass/density bin ranges from O.IZ© to below l.Oe- 
4 Z0. Again, ejection models have a greater effect on the 
distribution of metallicities than the reheating models, 
which is unsurprising, as the reheating portion of SAGE 
moves gas from the disk to the halo and the ejection portion 
moves it out to the IGM. 

The metallicities also show a second peak for very low 
mass galaxies in similar density environments. These galax¬ 
ies are below the resolution of the code, so this is most 
likely a numerical issue. Satellite galaxies do not have ejected 
reservoirs, as they are added to the host halo during infall, 
and therefore are not included in this figure. 

Focusing specifically on the mass and density range 
where the most metals are ejected, we define a subset of 
galaxies with 9 log^Q(M*/h“^ Mq) ^ 11 and 0.5 ^ 
logj^Q(p/p) ^ 2.0. By limiting the density, we eliminate the 
regions of space that are not dense enough to have mea¬ 
surable clouds of HI as well as higher density regions that 
preferentially host large groups and clusters where the hot 
halo dominates the gas content. Our stellar mass lower limit 



logio(p/p) 



Iogio(p/P) 

Figure 7. Average metallicity of the IGM around a galaxy as 
a function of stellar mass and environment. Top to bottom, the 
panels are RlEl, R2E2, and R3E3, respectively. The distributions 
are very similar for constant ejection models. 
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Figure 8. Distribution of metallicities in the IGM surround¬ 
ing galaxies in regions of density 0.5 Si logio(p/p) ^ 2.0 

and with stellar mass 9 Si logjg(M*/h“^ Mq) S 11. The 
dark/medium/light lines are reheating models 1, 2, and 3, re¬ 
spectively. The solid/dashed/dotted lines are ejection models 1, 
2, and 3, respectively. 


allows us to include only very well resolved galaxies. The 
high end cut is to remove galaxies larger than L*, which 
likely live in halos that are too massive to allow metals to 
escape. Galaxies in this sample account for 86 % of the met¬ 
als in the IGM in the fiducial model and 85-95% of galaxies 
in all of the other models. These numbers can be found in 
Table [T] 

Figure shows the distribution of metallicity measure¬ 
ments in the IGM surrounding galaxies in our subset. The 
models vary in both the value of their most likely metallic¬ 
ity (the peak of the distribution) as well as the number of 
low metallicity regions surrounding galaxies in the subset. 
The variation in scaling comes from the relative number of 
galaxies in this range, which can be seen in the SMFs in 
Figure 

The most likely metallicity (i.e. the peak in the distribu¬ 
tion) for each of the models is between —1.5 < logj^g(Z/Z 0 ) 
< —1.0. E2 models have slightly higher peak metallicities 
for a given reheating model. Despite ejecting fewer metals 
in general, the metals they do eject go into less dense regions 
(see Figure]^, decreasing the amount of IGM in the cell and 
increasing the metallicity. 

R1E2 and R2E2 have more low metallicity IGM regions 
than other models. This is partially due to the excess of 
galaxies produced by E2 models (middle panel. Figure]^ - 
a factor of 2 in this subset. These extra galaxies are below 
the mass threshold in the other models, and inhabit envi¬ 
ronments on the lower cusp of our density range. They eject 
small amounts of metal into lower density regions, doubling 
the number of regions in the IGM where metallicity can be 
measured, but having little effect (~ 10 %) on the relative 
amount of IGM metals found in the subset. So few of the 
galaxies in R3E2 eject metals that the increase in number of 
galaxies produced is entirely countered by the low percent¬ 
age of metal ejecting galaxies. 

The ejection models not only affect the amount of gas 


El 

E2 

E3 


% of Gals in Subset 


% of IGM Metals 


w/ Ejected Metals 


in Subset 


Rl 

R2 

R3 

Rl 

R2 

R3 

93% 

92% 

70% 

86% 

86% 

85% 

90% 

92% 

35% 

93% 

95% 

94% 

93% 

89% 

87% 

86% 

85% 

85% 


Table 1. Golumns 1-3: Percentage of galaxies with densities in 
the range 0.5 < logj^Q(p/p) < 2.0 and having a mass within 9.0 < 
Mq< 11.0 that have ejected metal reservoirs. Columns 
4-6: The percentage of total metals that are found in the IGM 
surrounding galaxies in this subset. 


El 

E2 

E3 


% of Metals Found 

Q 

z X 10-5 


in 

the IGM 

in 

the IGM 

Rl 

R2 

R3 

Rl 

R2 

R3 

8.1% 

9.5% 

4.6% 

0.835 

0.866 

0.383 

4.6% 

8.0% 

2.6% 

0.422 

0.390 

0.201 

8.1% 

7.4% 

6.8% 

0.835 

0.727 

0.420 


Table 2. Columns 1-3: The percentage of total metals that are 
found in the IGM. Columns 4-6: The total density of metals in 
the IGM. 


being ejected but also if gas is ejected at all. In models that 
use Rl, 90-93% of galaxies in this stellar mass and density 
range eject gas and metals, and in R2 models, the range is 
89-92%. In R3 models, however, the range is 35-87%; here, 
the presence of ejected gas is more tied to the reheating 
model than the ejection model. The percentages are listed in 
Table The relative numbers of galaxies that have ejected 
metals are apparent in Figure 

Since we do not change any parameters in SAGE other 
than the mass loading factors, we cannot rule out any of the 
models from just the stellar mass function and supernova 
feedback models. 


4 COMPARISONS TO OTHER WORK 
4.1 SWAMI vs. Hydrosimulations 


Theoretical models of the intergalactic medium tend to be 
relegated to hydrosimulations, so it is natural to begin by 
comparing our results to those. In Figure we show to¬ 
tal metallicity versus environment. By total metallicity, we 
mean the total mass of metals in a certain environment, di¬ 
vided by the total IGM gas mass in the same environment. 
It is comparable to Figure 4 in Oppenheimer et al. (2012 1 , 
who also find —1.5 < log]^Q(Z/Z 0 ) < —1.0 at p ~ 10 x p in 
all of their models with wind. 

We do not find the same high metallicity in high density 
environments that Oppenheimer et al. (20121 do for all of 
their models. This is mostly due to their inclusion of hot 
halo gas in their analysis. Larger galaxies reside in denser 
environments (Figure and in larger halos. Larger halos 
have more hot gas and a higher metallicity than the IGM in 
general because metals merely have to be pushed out of the 
disk rather than out of the halo entirely. 

The upturn we see at p > 300 x p is due to the lack of 
“pristine” IGM at high densities. Most of the dark matter 
particles in these environments are bound in halos, so the 
IGM has the metallicity of the ejected gas, rather than the 
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Figure 9. Total metallicity as a function of environments. The 
dark/medium/light lines are reheating models 1, 2, and 3, respec¬ 
tively. The solid/dashed/dotted lines are ejection models 1, 2, and 
3, respectively. 


metallicity of the mixed ejected-l-pristine gas found at lower 
densities. 

At low densities, where halos tend to be smaller and 
therefore the hot halo is less of a concern, we find the metal¬ 
licity drops off sharply. This is a good match to |Oppen-| 


heimer et al. ( 20121 ’s preferred model (vzw) of momentum 


conserved winds. Their constant wind (cw) model, which is 
similar to our R1 models, show a flat distribution in the 
metallicity in the range of p < p < 10 x p, and a sharp 
decrease at lower densities. 


Stinson et al. (20121 and Brook et al. (20121 also show 


metals outside the halo, with a sharp drop off at a few hun¬ 
dred kpc. Their 01V column density profiles indicate that 
the presence of metals is more of a smooth gradient through 
the halo into the IGM than our model is able to currently 
convey. We hope to explore this further in future work. 


4.2 SWAMI vs. Observations 


One of the many questions about metals in the Universe is 
whether they are all accounted for by current observations 
( [Pettini et al.|[l999P or is there a percentage that we are 
“missing”. This is a difficult question to answer due to un¬ 
certainty in the amount of metals that have been produced 
since the Big Bang ( |Madau &: Dickinson|20f^ and the dif¬ 
ficulty in measuring the quantity of metals in the IGM and 
hot halo gas. While theoretical models have the option of 
considering the IGM and its metals in any region of inter¬ 
est, current observational studies are limited to absorption 
lines. This makes one-to-one comparisons difficult, although 
the number of measurements is increasing with recent in¬ 
struments allowing for more statistical studies. 

We can use the SWAMI data to compare the net amount 
of metals in the IGM as well as the distribution of measured 
metallicities. Both observations ( Shull et al.||2014 l and hy¬ 
drosimulations ( Oppenheimer et al.|2012 | predict 10±5% of 
metals produced are in the IGM at z = 0. We hnd that in 5 
of our 8 unique models, M^(IGM)/M^ (Total) ~ 5-10%. Of 


the other three, two more (R1E2 and R3E1) are just below 
5%, leaving only R3E2 (where most of the supernova energy 
is used to heat the disk gas) with much too low a metal 
content. The values for each model can be found in Tabled 


4-2.1 Total metals in the IGM 


To address the question of the potentially missing metals. 


Danforth et al. (20141 have completed the largest IGM sur¬ 


vey to date, consisting of 75 AGN sight lines that pass 
through over 2,500 unique HI absorption systems. They find 
in a later paper ( Shull et al.||2014 | that the average metal¬ 
licity of the IGM in the low redshift universe is Q-z ~ 10“®. 

In comparison, we find Q-z = 0.2 — 0.9 x 10“®, depend¬ 
ing on the model, with E2 having, once again, the lowest 
metallicity for R1 and R3 reheating models. This is a direct 
result of limiting the supernova velocity in E2 to L* circu¬ 
lar velocity and below. E2 does not assume that all of the 
supernova energy is put into reheating and ejection, so it 
follows that it would eject less than its counterparts. 

Likewise, for a given ejection model, the cold gas- 
dependent reheating method (R3) gives the lowest IGM 
metal density. As discussed above in Section |3.2[ R3 has 
the highest disk mass loading factor for the mass range of 
interest. This means that reheating the gas takes more of 
the available energy and so less gas and fewer metals are 
ejected. 

All combinations of that include R1 or R2 are between 
fiz = 0.6 — 0.9 X 10“®, and are reasonable fits to the obser¬ 
vations. Moreover, our agreement with Shull et al.| ( |2014 rs 
observations is consistent with no “missing metals” problem. 
The values for each model can be found in Tabled 


4-2.2 Distribution of Metallicity 


The most commonly quoted value for the metallicity of the 
IGM is O.IZq, which comes from hydrosimulations such as 


Gen et al. (2001 1 . Measurements of quasar absorption lines 


range from —1.8 ^ [0/H] ^ —0.6 at z = 0.12 (Tripp et al. 
200l| to - 1.93 ^ [O/H] s; 0.03 at z = 0.08 - 0.23 ^Savage 
et al.||2014 ). Other measurements have similar ranges (e.g. 
Tumlinson et al.|2011| . 

Our maximum metallicities are all approximately so¬ 
lar, with R2E1 having the lowest maximum at logj^Q(Z/Z 0 ) 
~ —0.136. Since Savage et al. (20141 find an IGM metal¬ 


licity higher than this, it would appear this model is less 
likely than the others, although it does not completely rule 
it out. All models have their most likely metallicities in the 
range —1.5 < log]^Q(Z/Z 0 ) < —1 (Figure]^, with minimums 
at zero metallicity and maximums at solar. The observed 
metallicity ranges fall well within the models. 

[Lehner et al.| l |2013[ ) find a bimodal distribution in the 
metallicity in the circumgalactic medium around 28 Lyman 
limit systems. While they are surveying the gas within the 
halo and this paper addresses the gas outside, their results 
are nevertheless indicative that our mixing model ( 100 % 
mixing on relatively short timescales) is an oversimplifica¬ 
tion. The metallicities we find in Figure 8 have slightly lower 
maximum and significantly lower minimum values, but that 
is to be expected as gas flows out of the halo. 

Since observations are limited to quasar sightlines, mea¬ 
suring the IGM’s proximity to a galaxy of a given mass is 
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a rare event. Nevertheless, Kacprzak et al. (20141 find a 
M* = 1O^°'®M0 {h = 0.7) galaxy at z ~ 0.2 with a stellar 
metallicity of [O/H] = -0.21±0.08 and an IGM metallicity 
(58 kpc away) of [X/H] = -1.12±0.02. This separation most 
likely puts the measurement within the hot halo, but the 
authors determine the gas is infalling and therefore comes 
from the IGM. Additionally, Prochaska et al. (20111 find 
OIV out to 300 kpc around sub L* galaxies, well outside of 
the virial radii. While none of our models fail to reproduce 
these observations in IGM metallicities, runs which use R2 
tend to have higher galaxy metallicities, making them less 
likely to match Kacprzak et al. (20141’s galaxy overall. 

We would need to tune SAGE to each reheat¬ 
ing/ejection combination to absolutely rule out any model, 
but by comparing the current results to observation, we can 
identify weaknesses in the models. E2 is based on the logical 
assumption that not all supernova energy is converted into 
kinetic energy. Yet it under-predicts the median and total 
metals in the IGM. For E2 to successfully reproduce observa¬ 
tions, there needs to be an additional mechanism that pushes 
metals into the IGM around sub-L* galaxies. One possibility 
is that it preferentially pushes relatively high metallicity gas 
from the halo. 


5 SUMMARY 

In this paper, we have measured the effect of supernova feed¬ 
back on the contamination of the IGM by ejected metals. 
To do this, we have tested nine combinations of reheating 
and ejection models from supernova feedback. The three re¬ 
heating models assume different physics for the disk mass 
loading factor: no dependence on the host galaxy, inversely 
dependent on the circular velocity, and dependent on the 
distribution of cold gas in the disk. Likewise, the three ejec¬ 
tion models assume a constant halo mass loading factor with 
either a constant or a circular-velocity dependent wind, and 
a halo mass loading factor that is dependent on the amount 
of mass that is reheated. 

To model the IGM, we have developed a method 
for rapidly modelling baryons in the diffuse intergalactic 
medium akin to the semi-analytic models of galaxy forma¬ 
tion. It can be used in conjunction with the galaxy models to 
measure the effects of feedback on a galaxy’s environment. 
By varying the efficiency with which a supernova reheats gas 
from the disk and then ejects hot gas from the halo, we can 
make predictions concerning the metallicity of the IGM. 

We measure the metallicity of the IGM to be at its high¬ 
est around galaxies with 10.0 <logj^Q(M*//i“^ Mq)< 10.5 in 
environments about 10 x p. The actual value of this maxi¬ 
mum is highly dependent on the combination of reheating 
and ejection models used. Most (> 90%) of the metals in 
the IGM are ejected by galaxies with M*< Mq 

in environments with p < 30 x p, regardless of the reheat¬ 
ing/ejection model used. Our results indicate that the re¬ 
gions around small groups are ideal places to see metals in 
the IGM compared to both observations and hydrosimula¬ 
tions. 

In agreement with both observations and hydrosimula¬ 
tions, we find 5-10% of all metals are in the IGM. Our metal¬ 
licity ranges also compare favourably with observations of 
the IGM near individual galaxies. 


A key question in current IGM studies is the evolution 
of its metallicity with redshift. To answer this question, we 
need to fully couple our galaxy model, SAGE, with the dif¬ 
fuse gas model, SWAMI. This will allow us to leave ejected 
metals behind as galaxies move through the cells, and will 
be the subject of future work. 
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